The many-body localization phase transition 



O 

» 

o 
O 



c 

I 

CO 



c3 



i 

c 

o 
o 



> 

(N 
On 
OS 

d 
o 



Arijeet Pal 1 and David A. Huse 1 
'Physics Department, Princeton University, Princeton, NJ 08544, USA 

We use exact diagonalization to explore the many-body localization transition in a random-field 
spin-1/2 chain. We examine the correlations within each many-body eigenstate, looking at all 
high-energy states and thus effectively working at infinite temperature. For weak random field the 
eigenstates are thermal, as expected in this nonlocalized, "ergodic" phase. For strong random field 
the eigenstates are localized, with only short-range entanglement. We roughly locate the localization 
transition and examine some of its finite-size scaling, finding that this quantum phase transition at 
nonzero temperature might be showing infinite-randomness scaling with a dynamic critical exponent 
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I. INTRODUCTION 

As originally proposed in Anderson's seminal paper pj , 
an isolated quantum system of many interacting degrees 
of freedom with quenched disorder may be localized, and 
thus generically fail to approach local thermal equilib- 
rium, even in the limits of long time and large systems, 
and for energy densities well above the system's ground 
state. In the same paper, Anderson also treated the lo- 
calization of a single particle- like quantum degree of free- 
dom, and it is this single-particle localization, without 
interactions, that has received most of the attention in 
the half-century since then. Much more recently, Basko, 
et al. have presented a very thorough study of many- 
body localization with interactions at nonzero tempera- 
ture, and the topic is now receiving more attention; see 
e.g. [Ml. 

Many-body localization at nonzero temperature is a 
quantum phase transition that is of very fundamental 
interest to both many-body quantum physics and sta- 
tistical mechanics: it is a quantum "glass transition" 
where equilibrium quantum statistical mechanics breaks 
down. In the localized phase the system fails to ther- 
mally equilibrate. These fundamental questions about 
the dynamics of isolated quantum many-body systems 
are now relevant to experiments, since such systems can 
be produced and studied with strongly-interacting ultra- 
cold atoms 14]. And they may become relevant for cer- 
tain systems designed for quantum information process- 
ing |15| . Also, many-body localization may be underlying 
some highly nonlinear low-temperature current-voltage 
characteristics measured in certain thin films 11611 . 



II. THE MODEL 

Many-body localization appears to occur for a wide va- 
riety of particle, spin or q-bit models. Anderson's origi- 
nal proposal was for a spin system [l[ ; the specific simple 
model we study here is also a spin model, namely the 
Heisenberg spin-1/2 chain with random fields along the 
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where the static random fields hi are independent ran- 
dom variables at each site i, each with a probability dis- 
tribution that is uniform in [—h, h\. Except when stated 
otherwise, we take J = 1. The chains are of length 
L with periodic boundary conditions. This is one of 
the simpler models that shows a many-body localization 
transition. Since we will be studying the system's be- 
havior by exact diagonalization, working with this one- 
dimensional model that has only two states per site al- 
lows us to probe longer length scales than would be pos- 
sible for models on higher-dimensional lattices or with 
more states per site. We present evidence that at infinite 
temperature, /3 = l/T = 0, and in the thermodynamic 
limit, L — > oo, the many-body localization transition at 
h = h c = 3.5 ± 1.0 does occur in this model. The usual 
arguments that forbid phase transitions at nonzero tem- 
perature in one dimension do not apply here, since they 
rely on equilibrium statistical mechanics, which is ex- 
actly what is failing at the localization transition. We 
also present indications that this phase transition might 
be in an infinite-randomness universality class with an 
infinite dynamical critical exponent z — > oo. 

Our model has two global conservation laws: total en- 
ergy, which is conserved for any isolated quantum sys- 
tem with a time- independent Hamiltonian; and total S z . 
The latter conservation law is not essential for localiza- 
tion, and its presence may affect the universality class 
of the phase transition. For convenience, we restrict our 
attention to states with zero total S z . 

For simplicity, we consider infinite temperature, where 
all states are equally probable (and where the sign of the 
interaction J does not matter). The many-body local- 
ization transition also occurs at finite temperature; by 
working at infinite temperature we remove one parame- 
ter from the problem, and use all the eigenstates from the 
exact diagonalization (within the zero total S z sector) of 
each realization of our Hamiltonian. We see no reason 
to expect that the nature of the localization transition 
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differs between infinite and finite nonzero temperature, 
although it is certainly different at strictly zero temper- 
ature [O] ■ Note that this is a quantum phase transition 
that occurs at nonzero (even infinite) temperature. Like 
the more familiar ground-state quantum phase transi- 
tions, this transition is a sharp change in the properties 
of the many-body eigenstates of the Hamiltonian, as we 
discuss below. But unlike ground-state phase transitions, 
the many-body localization transition at nonzero temper- 
ature appears to be only a dynamical phase transition 
that is invisible in the equilibrium thermodynamics Q. 

There are many distinctions between the localized 
phase at large random field h > h c and the delocal- 
ized phase at h < h c . We call the latter the "ergodic" 
phase, although precisely how ergodic it is remains to 
be fully determined [l8j]. The distinctions between the 
two phases all are due to differences in the properties of 
the many-body eigenstates of the Hamiltonian, which of 
course enter in determining the dynamics of the isolated 
system. 

In the ergodic phase (h < h c ), the many-body eigen- 
states are thermal |19l - l22l ]. so the isolated quantum sys- 
tem can relax to thermal equilibrium under the dynam- 
ics due to its Hamiltonian. In the thermodynamic limit 
(L — > oo), the system thus successfully serves as its own 
heat bath in the ergodic phase. In a thermal eigen- 
state, the reduced density operator of a finite subsys- 
tem converges to the equilibrium thermal distribution 
for L — > oo. Thus the entanglement entropy between 
a finite subsystem and the remainder of the system is, 
for L — > oo, the thermal equilibrium entropy of the sub- 
system. At nonzero temperature, this entanglement en- 
tropy is extensive, proportional to the number of degrees 
of freedom in the subsystem. 

In the many-body localized phase (h > h c ), on the 
other hand, the many-body eigenstates are not thermal 
2): the "Eigenstate Thermalization Hypothesis" (l9M2^ 
is false in the localized phase. Thus in the localized 
phase, the isolated quantum system does not relax to 
thermal equilibrium under the dynamics of its Hamilto- 
nian. The infinite system fails to be a heat bath that can 
equilibrate itself. It is a "glass" whose local configura- 
tions at all times are set by the initial conditions. Here 
the eigenstates do not have extensive entanglement, mak- 
ing them accessible to DMRG-like numerical techniques 
[j| . A limit of the localized phase that is simple is J = 
with h > 0. Here the spins do not interact, all that hap- 
pens dynamically is local Larmor precession of the spins 
about their local random fields. No transport of energy 
or spin happens, and the many-body eigenstates are sim- 
ply product states with each spin either "up" or "down" . 

Any initial condition can be written as a density ma- 
trix in terms of the many-body eigenstates of the Hamil- 
tonian as p = J2mn Pmn\m){n\- The eigenstates have 
different energies, so as time progresses the off-diagonal 
density matrix elements m ^ n dephase from the par- 
ticular phase relations of the initial condition, while the 
diagonal elements p nn do not change. In the ergodic 



phase for L — > oo all the eigenstates are thermal so this 
dephasing brings any finite subsystem to thermal equi- 
librium. But in the localized phase the eigenstates are 
all locally different and athermal, so local information 
about the initial condition is also stored in the diagonal 
density matrix elements, and it is the permanence of this 
information that in general prevents the isolated quan- 
tum system from relaxing to thermal equilibrium in the 
localized phase. 

Our goals in this paper are (i) to present results in 
the ergodic and localized phases that are consistent with 
the expectations discussed above, and (ii), more impor- 
tantly, to examine some of the properties of the many- 
body eigenstates of our finite-size systems in the vicinity 
of the localization transition to try to learn about the na- 
ture of this phase transition. Although the many-body 
localization transition has been discussed by a few au- 
thors, there does not appear to be any proposals for the 
nature (the universality class) of this phase transition or 
for its finite-size scaling properties, other than some very 
recent initial ideas in Ref . [l(| • It is our purpose here to 
investigate these questions, extending the previous work 
of Oganesyan and Huse |4|, who looked at the many- 
body energy-level statistics of a related one-dimensional 
model. Since the many-body eigenstates have extensive 
entanglement on the ergodic side of the transition, it may 
be that exact diagonalization (or methods of similar com- 
putational "cost" [l(|) is the only numerical method that 
will be able to access the properties of the eigenstates on 
both sides of the transition. 



III. DOES IT THERMALIZE? 

As a first simple measure to probe how thermal the 
many-body eigenstates appear to be, we have looked at 
the local expectation value of the z component of the 
spin 

mtJ = {n\St\n) a (2) 

at site i in sample a in eigenstate n. For each site in 
each sample we compare this for eigenstates that are ad- 
jacent in energy, showing the mean value of the difference: 

[l m ic^ — m ia +1 ' > l] f° r various L and h in Fig. 1, where 
the eigenstates are labeled with n in order of their en- 
ergy. The square brackets denote an average over states, 
samples and sites. The number of samples used in the 
data shown in this paper ranges from 10 4 for L — 8, to 
50 for L — 16 and some values of h. In our figures we 
show one-standard-deviation error bars. Here and in all 
the data in this paper we restrict our attention to the 
many-body eigenstates that are in the middle one-third 
of the energy-ordered list of states for their sample. Thus 
we look only at high energy states and avoid states that 
represent low temperature. In this energy range, the dif- 
ference in energy density between adjacent states n and 
(n+1) is of order yf~L2~ L and thus exponentially small in 
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L as L is increased. If the eigenstates are thermal then 
adjacent eigenstates represent temperatures that differ 
only by this exponentially small amount, so the expecta- 
tion value of Sf should be the same in these two states 
for L — > oo. From Fig. 1, one can see that the differ- 
ences do indeed appear to be decreasing exponentially 
with increasing L in the ergodic phase at small h, as 
expected. [Here and throughout this paper, when we 
use logarithms, they are base e ("natural").] In the lo- 
calized phase at large h, on the other hand, the differ- 
ences between adjacent eigenstates remain large as L is 
increased, confirming that these many-body eigenstates 
are not thermal. 
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generic finite-size system has no degeneracies and the off- 
diagonal matrix elements of p each time-average to zero. 
As a result, the long-time average of the spin polarization 
in this mode is 



(Mi> 0o = |r^|M 1 t |ri)<n|M 1 |r 1 ) 



(5) 



Thus for each many-body eigenstate in each sample we 
can quantify how much it contributes to the initial and 
to the long-time averaged polarization. We then define 
the fraction of the contribution to the initial polarization 
that is dynamic and thus decays away (on average) at 
long time, as 



f, 



(n) 
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(6) 



In the ergodic phase, the system does thermalize, so the 
initial polarization does relax away and /<* 1 for L — > 
oo. In the localized phase, on the other hand, there is 
no long-distance spin transport, so f a n ' —> for L — >• oo. 
In Fig. 2 we show the mean values of / for each L vs. 
h. They show the expected behavior in the two phases 
(trending with increasing L towards either 1 or 0), and 
the phase transition is indicated by the crossover between 
large and small / that occurs more and more abruptly 
as L is increased. 



FIG. 1: (Color online) The logarithm of the mean difference 
between the local magnetizations in adjacent eigenstates (see 
text). The values of the random field h are indicated in the 
legend. In the ergodic phase (small h) where the eigenstates 
are thermal these differences vanish exponentially in L as L 
is increased, while they remain large in the localized phase 
(large h). 

Thcrmalization requires the transport of energy. In the 
present model with conserved total S z , it also requires 
the transport of spin. To study spin transport on the 
scale of the sample size L, we consider the relaxation of 
an initially inhomogeneous spin density: 

Mi =J2S-exp(i27rj/L) (3) 

3 

is the longest wavelength Fourier mode of the spin den- 
sity. Consider an initial condition that is at infinite 
temperature, but with a small modulation of the spin 
density in this mode, so the initial density matrix is 
po = (I + eM\)/Z, where e is infinitesimal, and Z is 
the partition function. The initial spin polarization of 
this mode is then 

(Mi)o = £>|/>oMi|n) = £ Y^{n\M\M x \n) . (4) 

n n 

If we consider a time average over long times, then 
the long-time averaged density matrix p^ is diagonal in 
the basis of the eigenstates of the Hamiltonian, since a 




FIG. 2: (Color online) The fraction of the initial spin po- 
larization that is dynamic (see text). The sample size L is 
indicated in the legend. In the ergodic phase (small h) the 
polarization decays substantially under the dynamics, while 
in the localized phase (large h) the decay is small, and this 
distinction gets sharper as L increases. 

A qualitatively similar finite-size scaling plot also in- 
dicating the phase transition is obtained by examin- 
ing the many-body eigenenergy spacings as was done 
in Ref. [3], and is shown as Fig. 3. We consider 

the level spacings 5 a n ^ = \E„ — si n_1 ^|, where 
is the many-body eigenenergy of eigenstate n in sam- 
ple a. Then we obtain the ratio of adjacent gaps as 
r a = mm\o a , da \/max{o a ,d a ), and average 
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this ratio over states and samples at each h and L. In 
the ergodic phase, the energy spectrum has GOE (Gaus- 
sian orthogonal ensemble) level statistics and the average 
value of r converges to [r] = 0.53 for L — > oo, while in 
the localized phase the level statistics are Poisson and 
[r] — >= 0.39. Note that our model is integrable at h = 0, 
so will not show GOE level statistics in that limit, and 
this effect is showing up for our smallest L and lowest h 
in Fig. 3. 
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FIG. 3: (Color online) The ratio of adjacent energy gaps (de- 
fined in the text). The sample size L is indicated in the legend. 
In the ergodic phase, the system has GOE level statistics, 
while in the localized phase the level statistics are Poisson. 

The crossings of the curves for different values of L 
in Figs. 2 and 3 give estimates of the location h c of 
the phase transition. Both plots show these estimates 
"drifting" towards larger h as L is increased, with the 
crossings at the largest L being slightly above h — 3. In 
both cases this "drifting" is also towards the localized 
phase, suggesting the behavior at the phase transition is, 
by these measures, more like the localized phase than it 
is like the ergodic phase. 

IV. SPATIAL CORRELATIONS 

To further explore the finite-size scaling properties of 
the many-body localization transition in our model, we 
next look at spin correlations on length scales of order 
the length L of our samples. One of the simplest correla- 
tion functions within a many-body eigenstate \n) of the 
Hamiltonian of sample a is 

C* z a (i,j) = (n\SfS?\n) a - (n\Sf\n) a (n\S?\n) a . (7) 

In Fig. 4 we show the mean value [log |C^(i, i + d)\] 
as a function of the distance d for representative values of 
h in the two phases and near the phase transition. Data 
are presented for various L. This correlation function 
behaves very differently in the two phases: 

In the ergodic phase, for large L this correlation func- 
tion should approach its thermal equilibrium value. For 
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FIG. 4: (Color online) The spin-spin correlations in the many- 
body eigenstates as a function of the distance d. The sample 
size L is indicated in the legend. The correlations decay ex- 
ponentially with d in the localized phase (h = 6.0), while they 
are independent of d at large d in the ergodic phase (h = 0.6). 
Intermediate behavior at h = 3.6, which is near the localiza- 
tion transition, is also shown. 



the states with zero total S z that we look at, (n\S z \n) = 
in the thermal eigenstates of the ergodic phase. However, 
the conservation of total S z does result in anticorrela- 
tions so that C^ z a {i,j) ~ — 1/(4(L— 1)) for well-separated 
spins. These distant spins at sites i and j are entangled 
and correlated: if spin i is flipped, that quantum of spin 
is delocalized and may instead be at any of the other 
sites, including the most distant one. These long-range 
correlations are apparent in Fig. 4 for h — 0.6, which 
is in the ergodic phase. Note that at large distance the 
correlations in the ergodic phase become essentially in- 
dependent of d — at large L and d, confirming that 
the spin flips are indeed delocalized. Although we only 
plot the absolute value of the correlations, in fact these 
correlations are almost all negative, as expected, in this 
large L ergodic regime. 

In the localized phase, on the other hand, the eigen- 
states are not thermal and (n\S z \n) remains nonzero for 
L — > oo. If spin i is flipped, within a single eigenstate 
that quantum of spin remains localized near site i, with 
its amplitude for being at site j falling off exponentially 
with the distance: C^(i,j) ~ exp (— \i — with £ 

the localization length. In the localized phase the typ- 
ical correlation and entanglement between two spins i 
and j thus fall off exponentially with the distance \i — j\ 
(except for \i ~ j\ near L/2, due to the periodic bound- 
ary conditions). This behavior is apparent in Fig. 4 for 
h = 6.0, which is in the localized phase and has a local- 
ization length that is less than one lattice spacing. We 
note that in the localized phase, as well as near the phase 
transition, the long distance spin correlations C zz are of 
apparently random sign. 

The data of Figs. 1-4 show the existence of and some of 
the differences between the ergodic and localiz ed p hases. 
We have also looked at entanglement spectra [23j of the 
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eigenstates (data not shown) , which also support the ro- 
bust existence of these two phases. In addition to con- 
firming the existence of these two distinct phases, we 
would like to locate and characterize the many-body lo- 
calization phase transition between them. However, in 
the absence of a theory of this transition, the nature of 
the finite-size scaling is uncertain, which makes it dif- 
ficult to draw any strong conclusions from these data 
with their modest range of L. In studies of ground-state 
quantum critical points with quenched randomness, very 
broadly speaking, one first step is to classify the transi- 
tions by whether they are governed (in a renormalization 
group treatment) by fixed points with finite or infinite 
randomness 0, 
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FIG. 6: (Color online) The scaled width a of the probability 
distribution of the logarithm of the long-distance spin corre- 
lations (see text). The legend indicates the sample lengths L. 
In the ergodic phase at small h and in the localized phase at 
large h, this width decreases with increasing L, while near the 
transition it increases. To produce the one-standard-deviation 
error bars shown, we have calculated the a (see text) for each 
sample by averaging only over sites and eigenstates within 
each sample, and then used the sample-to-sample variations 
of a to estimate the statistical errors. We have also (data not 
shown) calculated a by instead averaging (f> and (f> 2 over all 
samples; this produces scaling behavior for a that is qualita- 
tively the same as shown here, but with a somewhat larger 
in the localized phase and near the phase transition. 



FIG. 5: (Color online) The probability distributions of the 
logarithm of the long distance spin-spin correlation in the 
many-body eigenstates for sample size L = 16 and the values 
of the random field h is indicated in the legend. 

To explore this question for our system, we next look at 
the probability distributions of the long distance spin cor- 
relations. For quantum-critical ground states governed 
by infinite-randomness fixed points, these probability dis- 
tributions are found to be very broad [24]. In particular, 
we look at 

<t> = \og\C™ a {i,i+{L/2))\ , (8) 

whose probability distributions for L — 16 are displayed 
in Fig. 5 for various values of h. Note the distributions 
are narrow, as expected, in the ergodic phase and consis- 
tent with log-normal, as expected, in the localized phase. 
In between, in the vicinity of the apparent phase transi- 
tion, the distributions are quite broad and asymmetric. 

To construct a dimensionless measure of how these dis- 
tributions change shape as L is increased, we divide <p 
by its mean, defining rj — <j>/[<j)\. Then we quantify the 
width of the probability distribution of rj by the stan- 
dard deviation a = \Z[ri 2 ] — 1. This quantity is shown 
in Fig. 6 vs. h for the various values of L. By this 
measure, in both the ergodic and localized phases the 
distributions become narrower as L is increased, as can 
be seen in Fig. 6. This happens in the localized phase 



because although the mean of — grows linearly in L, 
the standard deviation is expected to grow only ~ \[L. 
Over the small range of L that we can explore, a is found 
to decrease more slowly than the expected L^ 1 / 2 in the 
localized phase, but it does indeed decrease. 

This scaled width cr^Qi) of the probability distribu- 
tion of as a function of the random field h for each 
sample size L shows a maximum between the ergodic 
and localized phases. In the vicinity of the phase tran- 
sition, a actually increases as L is increased, suggesting 
that its critical value is nonzero, like for quantum-critical 
ground states that are governed by an infinite random- 
ness fixed point. This suggests the possibility that this 
one-dimensional many-body localization transition might 
also be in an infinite-randomness universality class. The 
peak in this plot is close to h = 4, and is thus suggesting 
a slightly higher estimate of h c than the crossings in Figs. 
2 and 3. 



V. DYNAMICS 

In the study of the spectral and localization proper- 
ties of noninteracting particles in finite samples (such 
as quantum dots), there are two very important energy 
scales: the level spacing 5 and the Thouless energy E^p . 
The Thouless energy is h times the rate of diffusive re- 
laxation on the scale of the sample. The diffusive (non- 
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localized or ergodic) phase is where Et is larger than 
6, and for d-dimensional samples with d > 3, the local- 
ization transition occurs when these two energy scales 
are comparable. Since the single-particle level spacing 
in a d-dimcnsional system of linear size L behaves as 
5 ~ L~ d and this sets the relaxation time at the local- 
ization transition, the dynamic critical exponent for the 
single-particle localization transition is z = d. 

A possibility that we will now investigate is that the 
many-body localization transition also occurs when the 
Thouless energy is of order the many-body level spac- 
ing. Since the many-body level spacing behaves as 
log 6 ~ — L d , this corresponds to an infinite dynamic crit- 
ical exponent z — > oo. Note also that even for our model 
with d = 1 this is a stronger divergence of the critical 
time scales than occurs at the known infinite-randomness 
ground-state quantum critical points, where log 6 ~ — lr 
with V < 1/2. 

0.25 1 1 1 , 1 
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FIG. 7: (Color online) Contribution to the dynamic part of 
(Mi) from matrix elements between adjacent energy states 
(see text). In the ergodic and localized phase the contribution 
is decreasing to zero with increasing sample size. The sample 
size L is indicated in the legend. The maximum contribution 
from adjacent states is close to the critical point. 

It is important to note that the model (1) we study has 
two globally conserved quantities; total energy and total 
S z . Their respective transport times (and hence their 
corresponding Thouless energy) in the ergodic phase may 
have different scaling properties close to the critical point. 
By studying the relaxation of the spin modulation, Mi , 
we are specifically probing the spin transport time which 
may diverge differently from the energy transport time 
close to the critical point. Such a possibility has been dis- 
cussed in the context of zero-temperature metal-insulator 
transitions (26j and may play a role in deciding the uni- 
versality class of the many-body localization transition. 

Naively, the Thouless energy is set by the relaxation 
rate of the longest-wavelength spin density modulation, 
Mi. If the scaling at the many-body localization tran- 
sition is such that the Thouless energy is of order the 
many-body level spacing, then at the transition a nonzero 
fraction of the dynamic part of (Mi) should be from its 



matrix elements between adjacent energy levels, and this 
fraction should remain large as L is increased. In each 
sample a, the contribution of a given eigenstate \n) to 
the dynamic part of (Mi) is given by 

(AMi)W = (n|MjMi|n) - |(n|Mi|n)| 2 . (9) 

In the ergodic phase, (AMi)i™' 1 has significant contribu- 
tions from matrix elements with many other eigenstates, 
and the Thouless energy is a measure of the energy range 
over which these contributions occur. To quantify this, 
we define as the contribution to the dynamic part 
of (Mi) from the matrix elements between state n and 
states n ± i: 

Qia = \(n-i\M 1 \n)\ 2 + \(n\M 1 \n + i)\ 2 (10) 

in sample a. Note that 

S^oQl: 3 = (AMi)W . (II) 

We define P^ n) = Q { "J/(AM 1 )t ) as the fraction of the 
longest-wavelength "diffusive" dynamics that is due to 
interference between adjacent (i = 1) many-body energy 
levels. Fig. 7 shows this quantity averaged over disorder 
realizations and states. 

If at the localization transition the Thouless energy Et 
is proportional to the many-body level spacing 5, then [P] 
should remain nonzero in the limit L — > oo . We do indeed 
find a strong peak in this fraction near the many-body 
localization transition, and that the fraction is large and 
not decreasing much as L is increased. Note that the 
level spacing decreases by almost a factor of 4 for every 
increase of L by two additional spins, so near the tran- 
sition the Thouless energy is apparently decreasing by 
almost the same factor as L is increased. This seems at 
least consistent with Et ~ 8 scaling, and thus dynamic 
exponent z — > oo. In the localized phase, the dynamics 
is due to spin-moves that are short-range in real space 
(probably of order the localization length). These spin- 
hops involve pairs of many-body eigenstates that become 
far apart (large i) for large L; this is why [P] drops with 
increasing L in the localized phase. Note that the peak 
in [P] occurs a little below h = 3. If one ignores L = 8, 
the location of this peak is apparently drifting to larger 
h with increasing L, consistent with our other rough es- 
timates of h c . 

The dynamic fraction [fa] (Fig. 2) tends to 1 in the 
ergodic phase and decreases to in the localized phase. 
The probability distribution of /!"' (P(f)) is strongly 
peaked around 1 and in these respective phases. At the 
phase transition, this distribution could either be peaked 
at the critical f c , broadly distributed, or even bimodal 
with peaks near both zero and one. In Fig. 8, we show 
P(f) for a disorder strength h = 3.0 close to the esti- 
mated transition, for system sizes 10 and 16. This dis- 
tribution P{f) becomes broader and more bimodal with 
increasing L. This feature of the distribution is consis- 
tent with the indication from Fig. 6 that the critical 
point may be governed by a strong disorder fixed point. 



7 




' 0.2 0.4 0.6 0.8 1 



f 

FIG. 8: (Color online) Probability distribution of the dynamic 
fraction of (Mi) for L — 10 and 16. Close to the transition for 
h — 3.0, the distribution becomes broader and more bimodal 
with increasing L. 

VI. SUMMARY 

This study of the exact many-body eigenstates of our 
model (1) has demonstrated some of the properties of the 
ergodic and localized phases. We also find a rough esti- 
mate of the localization transition using various different 



diagnostics. Based on earlier work by one of the authors 
[Ij , the many-body energies go from having GOE to Pois- 
son level statistics with increasing disorder. The scaling 
of the probability distributions of the long-distance spin 
correlations suggests that the transition might be gov- 
erned by an infinite-randomness fixed point with dynamic 
critical exponent z —¥ oo. We also study the relaxation of 
spin modulation under the dynamics of the Hamiltonian. 
In this case our results are consistent with Et ~ S scaling 
at criticality, in apparent agreement with our earlier con- 
clusion of z — » oo at the transition. These results suggest 
that efforts to develop a theory of this interesting phase 
transition should consider the possibility of a strong dis- 
order renormalization group approach. Of course, the 
model we have studied is only one-dimensional, and the 
behavior of this transition in higher dimensions might be 
different in important ways. 
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